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Data from the LIGO Livingston interferometer and the ALLEGRO resonant bar detector, taken 
during LIGO's fourth science run, were examined for cross-correlations indicative of a stochastic 
gravitational-wave background in the frequency range 850-950 Hz, with most of the sensitivity arising 
between 905 Hz and 925 Hz. ALLEGRO was operated in three different orientations during the 
experiment to modulate the relative sign of gravitational-wave and environmental correlations. No 
statistically significant correlations were seen in any of the orientations, and the results were used 
to set a Bayesian 90% confidence level upper limit of f2 gw (/) < 1.02, which corresponds to a 
gravitational wave strain at 915 Hz of 1.5 x 10~ 23 Hz -1 / 2 . In the traditional units of /iioo^gw(/), this 
is a limit of 0.53, two orders of magnitude better than the previous direct limit at these frequencies. 
The method was also validated with successful extraction of simulated signals injected in hardware 
and software. 

PACS numbers: 



I. INTRODUCTION 

One of the signals targeted by the current genera- 
tion of ground-based gravitational wave (GW) detectors 
is a stochastic gravitational-wave background (SGWB) 
[2 9 • Such a background is analogous to the cosmic 
microwave background, although the dominant contribu- 
tion is unlikely to have a blackbody spectrum. A SGWB 
can be characterized as cosmological or astrophysical in 
origin. Cosmological backgrounds can arise from, for 
example, pre-big-bang models [1, [H, [|| , amplification of 
quantum vacuum fluctuations during inflation [7|, [a IS | , 




phase transitions [10L llll ] , and cosmic strings 
Astrophysical backgrounds consist of a superposition of 
unresolved sources, which can include rotating neutron 
star s [TBI . [n| , supernovae [TtJ and low-mass X-ray bina- 
ries |l8[. 

The standard cross-correlation search for a SGWB 
necessarily requires two or more GW detectors. Such 
searches have been performed using two resonant bar de- 
tectors [2(| and also using two or more kilometer-scale 
GW interferometers (IFOs) [H, [H, [H]. The present 
work describes the results of the first cross-correlation 
analysis carried out between an IFO (the 4 km IFO at 
the LIGO Livingston Observatory (LLO), known as LI) 
and a bar (the cryogenic ALLEGRO detector, referred to 
as Al). This pair of detectors is separated by only 40 km, 
the closest pair among modern ground-based GW detec- 
tor sites, which allows it to probe the stochastic GW 
spectrum around 900 Hz. In addition, the ALLEGRO 
bar can be rotated, changing the response of the corre- 
lated data streams to stochastic GWs and thus providing 
a means to distinguish correlations due to a SGWB from 
those due to correlated environmental noise [25[ . This pa- 
per describes cross-correlation analysis of LI and Al data 
taken between February 22 and March 23, 2005, during 
LIGO's fourth science run (S4). Average sensitivities of 



a Member of ALLEGRO Collaboration 

b Member of LIGO Scientific Collaboration and ALLEGRO Collab- 
oration 



LI and Al during S4 are shown in Fig. [T] ALLEGRO 
was operated in three orientations, which modulated the 
GW response of the LLO- ALLEGRO pair through 180° 
of phase. 
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FIG. 1: Sensitivity of the LLO IFO (LI) and the AL- 
LEGRO bar (Al) during S4, along with strain associated 
with f2 gw (/) = 1.02 (assuming a Hubble constant of Ho = 
72 km/s/Mpc). (There are two fi gw (/) = L02 curves, cor- 
responding to the different strain levels such a background 
would generate in an IFO and a bar, as explained in SeclTTl 
and [27J •) The quantity plotted is amplitude spectral den- 
sity (ASD), the square root of the one-sided power spectral 
density defined in (|4.2|l . at a resolution of 0.25 Hz. 



The LLO-ALLEGRO correlation experiment is com- 
plementary to experiments using data from the two LIGO 
sites, in that it is sensitive to a SGWB at frequencies 
of around 900 Hz rather than 100 Hz. Targeted sources 
are thus those with a relatively narrow-band spectrum 
peaked near 900 Hz. Spectra with such shapes can arise 
from exotic cosmological models, as described in Sec. [Til 
or from astrophysical populations (l6| . 

The organization of this paper is as follows. Sec. HT1 
reviews the properties and characterization of a SGWB. 
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Sec. IIIII describes the LLO and ALLEGRO experimen- 
tal arrangements, including the data acquisition and 
strain calibration for each instrument. Sec. IIVI describes 
the cross-correlation method and its application to the 
present situation. Sec.[V]describes the details of the post- 
processing methods and statistical interpretation of the 
cross-correlation results. Sec. I VII describes the results of 
the cross-correlation measurement and the correspond- 
ing upper limit on the SGWB strength in the range 850- 
950 Hz. Sec. IVIII describes the results of our analysis 
pipeline when applied to simulated signals injected both 
within the analysis software and in the hardware of the 
instruments themselves. Sec. IVIIII compares our results 
to those of previous experiments and to the sensitivities 
of other operating detector pairs. Sec. IIXI considers the 
prospects for future work. 



II. STOCHASTIC GRAVITATIONAL- WAVE 
BACKGROUNDS 

A gravitational wave (GW) is described by the metric 
tensor perturbation h a b(f^t). A given GW detector, lo- 
cated at position r^et on the Earth, will measure a GW 
strain which, in the long-wavelength limit, is some pro- 
jection of this tensor: 



Overlap Reduction Function 



h(t) = hab(fdet,t)d 



a b 



(2.1) 



where d ab is the detector response tensor, which is 

1 



d$ o) = -(x a x b -y a y b ) 



(2.2) 



for an interferometer with arms parallel to the unit vec- 
tors x and y and 



"(ba 



* a *b 
U U 



(2.3) 



for a resonant bar with long axis parallel to the unit 
vector u. 

A stochastic GW background (SGWB) can arise from a 
superposition of uncorrelated cosmological or astrophys- 
ical sources. Such a background, which we assume to be 
isotropic, unpolarizcd, stationary and Gaussian, will gen- 
erate a cross-correlation between the strains measured by 
two detectors. In terms of the continuous Fourier trans- 
form defined by a(f) = dta{t) exp(— i2Trft), the ex- 
pected cross-correlation is 

(K(f)Mf')) = \s(f-f)S gw (f) 112(f) (2.4) 
where 
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FIG. 2: The overlap reduction function for LIGO Livingston 
Observatory (LLO) with ALLEGRO and with LIGO Han- 
ford Observatory (LHO). The three LLO-ALLEGRO curves 
correspond to the three orientations in which ALLEGRO 
was operated during LIGO's S4 run: "XARM" (N108°W) is 
nearly parallel to the x-arm of LLO ("aligned"); "YARM" 
(N18°W) is nearly parallel to the y-arm of LLO ("anti- 
aligned"); "NULL" (N63°W) is halfway in between these 
two orientations (a "null alignment" midway between the two 
LLO arms). Note that for non-zero frequencies, the separa- 
tion vector between the two sites breaks the symmetry be- 
tween the "XARM" and "YARM" alignments, and leads to 
an offset of the "NULL" curve, as described in JQ. The LLO- 
LHO overlap reduction function is shown for reference. The 
frequency band of the present analysis, 850 Hz < / < 950 Hz, 
is indicated with dashed vertical lines. 



vector h. The ORF for several detector pairs of interest 
is shown in Fig. 

S g vj(f) is the one-sided spectrum of the SGWB. This 
is the one-sided power spectral density (PSD) the back- 
ground would generate in an interferometer with per- 
pendicular arms, which can be seen from (|2.4|) and the 
fact that the ORF of such an interferometer with itself 
is unity. Since the ORF of a resonant bar with itself is 
4/3 (see [27| and Sec.EH]for more details), the PSD of 
the strain measured by a bar detector due to the SGWB 
would be (4/3)5gw(/). 

A related measure of the spectrum is the dimensionlcss 
quantity f2 gw (/), the GW energy density per unit loga- 
rithmic frequency divided by the critical energy density 
p c needed to close the universe: 



7i 2 (/) = d lab d c 2 d ^JJ d*n A e «*/A-(iWi)/c 

n (2 ' 5) 

is the overlap reduction function (ORF) [26| between the 
two detectors, defined in terms of the projector -P TTn ^ 
onto traceless symmetric tensors transverse to the unit 



O < f \ — — 

Pc df 



IOtt 2 
3H$ 



f 3 S gw (f) 



(2.6) 



Note that the definition £l gw (/) thus depends on the 
value of the Hubble constant Hq. Most SGWB literature 
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avoids this artificial uncertainty by working in terms of 

^ = ( ioow,/ M pc ) ,n ^ (27) 

rather than O gw (/) itself. We will instead follow the 
precedent set by |23j and quote numerical values for 
O gw (/) assuming a Hubble constant of 72km/s/Mpc. 

A variety of spectral shapes have been proposed for 
f2 gw , for both astrophysical and cosmological stochastic 
backgrounds 0, [3(|[3l[. For example, whereas the slow- 
roll inflationary model predicts a constant Q gw (/) in the 
bands of LIGO or ALLEGRO, certain alternative cosmo- 
logical models predict broken-power law spectra, where 
the rising and falling slopes, and the peak- frequency are 
determined by model parameters 0] . String- inspired pre- 
big-bang cosmological models belong to this category 
[S HI]- F° r certain ranges of these three parameters 
the LLO- ALLEGRO correlation measurement offers the 
best constraints on theory that can be inferred from any 
contemporary observation. This can happen, e.g., if the 
power-law exponent on the rising spectral slope is greater 
than 3 and the peak-frequency is sufficiently close to 900 

Hz m. 

III. EXPERIMENTAL SETUP 
A. The LIGO Livingston Interferometer 

The experimental setup of the LIGO observatories has 
been described at length elsewhere 38]. Here we pro- 
vide a brief review, with particular attention paid to de- 
tails significant for the LLO- ALLEGRO cross-correlation 
measurement. 

The LIGO Livingston Observatory (LLO) is an inter- 
feromctric GW detector with perpendicular 4-km arms. 
The laser interferometer senses directly any changes in 
the differential arm length. It does this by splitting a 
light beam at the vertex, sending the separate beams 
into 4-km long optical cavities of their respective arms, 
and then recombining the beams to detect any change in 
the optical phase difference between the arms, which is 
equivalent to a difference in light travel time. This pro- 
vides a measurement of h(t) as defined in (|2.ip and (|2.2j) . 
However, the measured quantity is not exactly h(t) for 
two reasons: 

First, there are local forces which perturb the test 
masses, and so produce changes in arm length. There 
are also optical and electronic fluctuations that mimic 
real strains. The combination of these effects causes a 
strain noise n(t) to always be present in the output, pro- 
ducing a measurement of 

s(t) = h(t) + n{t) (3.1) 

Second, the test masses are not really free. There is 
a servo system, which uses changes in the differential 
arm length as its error signal q(t), and then applies extra 



("control") forces to the test masses to keep the differ- 
ential arm length nearly zero. It is this error signal q(t) 
which is recorded, and its relationship to the strain esti- 
mate s(t) is most easily described in the Fourier domain: 

«(/) - R(fMf) (3-2) 

The response function R(f) is estimated by a combina- 
tion of modelling and measurement (39l | and varies slowly 
over the course of the experiment. 

Because the error signal q(t) has a smaller dynamic 
range than the reconstructed strain s(t), our analysis 
starts from the digitized time series q[k] = q(tk) (sam- 
pled 2 14 = 16384 times per second, and digitally down- 
sampled to 4096 Hz in the analysis) and reconstructs the 
LLO strain only in the frequency domain. 

B. The ALLEGRO Resonant Bar Detector 

The ALLEGRO resonant detector, op erated by a 
group from Louisiana State University [4l|, is a two- ton 
aluminum cylinder coupled to a niobium secondary res- 
onator. The secondary resonator is part of an inductive 
transducer [42| which is coupled to a DC SQUID. Strain 
along the cylindrical axis excites the first longitudinal 
vibrational mode of the bar. The transducer is tuned 
for sensitivity to this mechanical mode. Raw data ac- 
quired from the detector thus reflect the high Q resonant 
mechanical response of the system. A major technical 
challenge of this analysis is due to the extent to which 
the bar data differ from those of the interferometer. 



1. Data Acquisition, Heterodyning and Sampling 

The ALLEGRO detector has a relatively narrow sen- 
sitive band of ~ 100 Hz centered around ~ 900 Hz near 
the two normal modes of the mechanical bar-resonator 
system. For this reason, the output of the detector can 
first be heterodyned with a commercial lock-in ampli- 
fier to greatly reduce the sampling rate, which is set at 
250 samples/s. Both the in-phase and quadrature out- 
puts of the lock-in are recorded and the detector output 
can thus be represented as a complex time series which 
covers a 250 Hz band centered on the lock-in reference 
oscillator frequency. This reference frequency is chosen 
to be near the center of the sensitive band, and during 
the S4 run it was set to 904 Hz. The overall timing of 
data heterodyned in this fashion is provided by both the 
sampling clock and the reference oscillator. Both time 
bases were locked to GPS. 

The nature of the resonant detector and its data acqui- 
sition system gives rise to a number of timing issues: het- 
erodyning, hltcr delays of the electronics and the timing 
of the data acquisition system itself [2l[ . It is of critical 
importance that the timing be fully understood so that 
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the phase of any potential signal may be recovered. Con- 
vincing evidence that all of the issues are accounted for 
in is demonstrated by the recovery and cross-correlation 
of test signals simultaneously injected into both detec- 
tors. The signals were recovered at the expected phase 
as presented in Sec. IVIII 



2. Strain Calibration 

The raw detector output is proportional to the dis- 
placement of the secondary resonator, and thus has a 
spectrum with sharp line features due to the high-Q res- 
onances of the bar-resonator system as can be seen in 
Fig. [3l The desired GW signal is the effective strain on 



Raw spectrum, S4 run - averaged over orientation 
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FIG. 3: The graph displays the amplitude spectral density 
of raw ALLEGRO detector output during S4, at a frequency 
resolution of 0.1 Hz. For this graph these data have not been 
transformed to strain via the calibrated response function The 
vertical scale represents digital counts/\/Hz The normal me- 
chanical modes where the detector is most sensitive are at 
880.78 Hz and 917.81 Hz. There is an injected calibration line 
at 837 Hz. Also prominent are an extra mechanical resonance 
at 885.8 Hz and a peak at 904 Hz (DC in the heterodyned data 
stream) 

the bar, and recovering this means undoing the resonant 
response of the detector. This response has a long co- 
herence time - thus long stretches of data are needed 
to resolve the narrow lines in the raw data. The strain 
data have a much flatter spectrum, as shown in Fig. [1] 
Therefore it is practical to generate the calibrated strain 
time series, s(t), and use that as the input to the cross 
correlation analysis. 

The calibration procedure, described in detail in [2lj . 
is carried out in the frequency domain and consists of 
the following: A 30 minute stretch of clean ALLEGRO 
data is windowed and Fourier-transformed. The mechan- 
ical mode frequencies drift slightly due to small temper- 



ature variations, so these frequencies are determined for 
each stretch and those are incorporated into the model of 
the mechanical response of the system to a strain. The 
model consists of two double poles at these normal mode 
frequencies. In addition to this response we must then 
account for the phase shifts due to the time delays in the 
lock-in and anti-aliasing niters. 

After applying the full response function, the data are 
then inverse Fourier-transformed back to the time do- 
main. The next 50% overlapping 30 minute segment is 
then taken. The windowed segments are stitched to- 
gether until the entire continuous stretch of good data 
is completed. The first and last 15 minutes are dropped. 
The result represents a heterodyned complex time series 
of strain, whose amplitude spectral density is shown in 

Fig.rrj 

The overall scale of the detector output in terms of 
strain is determined by applying a known signal to the 
bar. A force applied to one end of the bar has a sim- 
ple theoretical relationship to an equivalent gravitational 
strain [2l|, [32|, [33[ . A calibrated force can be applied via 
a capacitive "force generator" which also provides the 
mechanism used for hardware signal injections. A recip- 
rocal measurement - excitation followed by measurement 
with the same transducer - along with known properties 
of the mechanical system, allows the determination of 
the force generator constant . With that constant deter- 
mined (with units of newtons per volt) a calibrated force 
is applied to the bar and the overall scale of the response 
determined. 



3. Orientation 

A unique feature of this experiment is the ability to ro- 
tate the ALLEGRO detector and modulate the response 
of the ALLEGRO-LLO pair to a GW background 
Data were taken in three different orientations of ALLE- 
GRO, known as XARM, YARM, and NULL, detailed in 
Table [U As shown in Fig. [5] and (|2.4p . these orientations 
correspond to different pair responses due to different 
overlap reduction functions. In the XARM orientation- 
the bar axis parallel to the x-arm of the interferometer- 
a GW signal produces positive correlation between the 
data in the two detectors detectors. In the YARM orien- 
tation a GW signal produces an anti-correlation. In the 
NULL oricntation-the bar aligned halfway between the 
two arms of the interferometer-the pair has very nearly 
zero sensitivity as a GW signal produces almost zero cor- 
relation between the detectors' data. A real signal is thus 
modulated whereas many types of instrumental correla- 
tion would not have the same dependence on orientation. 



IV. CROSS-CORRELATION METHOD 

This section describes the method to used to search 
for a SGWB by cross-correlating detector outputs. In the 
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Dates 


orientation 


azimuth 


7(850 Hz) 


7(915 Hz) 


7(950 Hz) 


2005 Feb 22-2005 Mar 4 


YARM 


N108°W 


-0.9087 


-0.8947 


-0.8867 


2005 Mar 4-2005 Mar 18 


XARM 


N18°W 


0.9596 


0.9533 


0.9498 


2005 Mar 18-2005 Mar 23 


NULL 


N63°W 


0.0280 


0.0318 


0.0340 



TABLE I: Orientations of ALLEGRO during the LIGO S4 Run, including overlap reduction function evaluated at the extremes 
of the analyzed frequency range, and at the frequency of peak sensitivity. Note that while the NULL orientation represents 
perfect misalignment (7 = 0) at Hz, it is not quite perfect at the frequencies of interest. This is primarily because of an 
azimuth- independent offset term in j(f) which contributes at non-zero frequencies [25l. |27([. Due to this term, it is impossible 
to orient ALLEGRO so that 7(/) = at all frequencies, and to set it to zero around 915 Hz one would have to use an azimuth 
of N62°W rather than N63°W. 



case of Ll-Al correlation measurements, it is complicated 
by the different sampling rates for the two instruments 
and the fact that the Al data are heterodyned at 904 Hz 
prior to digitization. 



A. Continuous-Time Idealization 

Both ground-based interferometric and resonant-mass 
detectors produce a time-series output which can be re- 
lated to a discrete sampling of the signal 



Si {t) = hi(t) + ni(t) 



(4.1) 



where i labels the detector (1 or 2 in this case), hi(t) is 
the gravitational- wave strain defined in (|2.1j) , and rij(i) 
is the instrumental noise associated with each detector, 
converted into an equivalent strain. The detector output 
is characterized by its power spectral density Pi(f) 



{^(f)Mf')) = ^(f-f , )p i (f) 



(4.2) 



which should be dominated by the auto-correlation of the 
noise ((s£ (/)£(/')) « (n*(f)ni (/')})• If the instrument 
noise is approximately uncorrelated, the expected cross- 
correlation of the detector outputs is [cf. (|2.4[) ] 



(4.3) 



= 2*(/-/')<V(/)7i2(/) 

which can be used along with the auto-correlation (|4.2p 
to determine the statistical properties of the cross- 
correlation statistic defined below. 

We use the optimally-filtered cross-correlation method 
described in [3 [22j to calculate a cross-correlation statis- 
tic which is an approximation to the continuous-time 
cross-correlation statistic 



Y c = / dt 1 dt 2 s 1 (t 1 )Q(t 1 -t 2 )s 2 (t 2 ) 



dfs\{f)Q{f)s 2 {f) 



(4.4) 



In the continuous-time idealization, such a cross- 
correlation statistic, calculated over a time T, has an 



expected mean 

HY- = (Y c ) « \ 
and variance 



df 7(1/1) S gw (f)Q(f) (4.5) 



/J.Y" 



T 
1 



dfPi(f)P 2 (f) Q(f) 



(4.6) 

Using (|4.5|) and (|4.6j) . the optimal choice for the filter 
Q(/), given a predicted shape for the spectrum S l gw (/) 
can be shown fijl to be 



Q(/)oc 



7(1/1) ggw(/) 



(4.7) 



If the spectrum of gravitational waves is assumed to be 
a power law over the frequency band of interest, a con- 
venient parameterization of the spectrum, in terms of 
Q gw (f) defined in 



(4.8) 



^gw(Z) — 7- 
, JR 



where fn us a conveniently-chosen reference frequency 
and £Ir = il gw (/fl)- The cross-correlation measurement 
is then a measurement of and if the optimal filter is 
normalized according to 



where 



jV = 



20n 2 



lif) if/far 

|/| 3 iM/)W) 

df bu) (f/f R r? 

P Pl(f)P2(f) 



(4.9a) 



(4.9b) 



then the expected statistics of Y c in the presence of a 
background of actual strength flu are 



/iyc =Q, R T 



(4.10) 



and 



T 



10tt 



2 \ z 



df w) {f/f R r\ 2 



,3£T V \J-00P A(/)W) 

(4.11) 

and a measurement of Y c /T provides a point estimate of 
the background strength flu with associated estimated 
errorbar ay? /T. 
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B. Discrete-Time Method 

1. Handling of Different Sampling Rates and Heterodyning 

Stochastic-background measurements using pairs of 
LIGO interferometers [12] have implemented (|4.4[) from 
discrete samplings s jfc] = s(t +k St) as follows: First the 
continuous Fourier transforms s(/) were approximated 
using the discrete Fourier transforms of windowed and 
zero-padded versions of the discrete time series; then 
an optimal filter was constructed using an approxima- 
tion to (|4.7j) . and finally the product of the three was 
summed bin-by-bin to approximate the integral over fre- 
quencies. The discrete version of Q(f) was simplified in 
two ways: first, because of the averaging used in calculat- 
ing the power spectrum, the frequency resolution on the 
optimal filter was generally coarser than that associated 
with the discrete Fourier transforms of the data streams, 
and second, the value of the optimal filter was arbitrar- 
ily set to zero outside some desired range of frequencies 
imin < / < /max- This was justified because the optimal 
filter tended to have little support for frequencies outside 
that range. 

The present experiment has two additional compli- 
cations associated with the discretization of the time- 
series data. First, the Al data are not a simple time- 
sampling of the gravitational-wave strain, but are base- 
banded with a heterodyning frequency f 2 = 904 Hz as 
described in Sec. IIIIB lllnTB2l Second, the Al data 
are sampled at (6*2 ) _1 = 250 Hz, while the LI data are 
sampled at 16384 Hz, and subsequently downsamplcd to 
(Sti)^ 1 = 4096Hz. This would make a time-domain 
cross-correlation extremely problematic, as it would ne- 
cessitate a large variety of time offsets t\ — t 2 . In the 
frequency domain, it means that the downsampled LI 
data, once calibrated, represent frequencies ranging from 
-2048 Hz to 2048 Hz, while the calibrated Al data rep- 
resent frequencies ranging from (904 — 125) Hz = 779 Hz 
to (904+ 125) Hz = 1029 Hz. These different frequency 
ranges do not pose a problem, as long as the range of fre- 
quencies chosen for the integral satisfies / m i n > 779 Hz 
and / m ax < 1029 Hz. Another requirement is that for the 
chosen frequency resolution, the Al data heterodyne ref- 
erence frequency must align with a frequency bin in the 
LI data. This is satisfied for integer-second data spans 
and integer-hertz reference frequencies. With these con- 
ditions, the Fourier transforms of the Al and LI data 
are both defined over a common set of frequencies, as 
detailed in [36[ . Looking at the Al spectrum in Fig. [TJ a 
reasonable range of frequencies should be a subset of the 
range 850 Hz < / < 950 Hz. 



• For LI, a real time series = 0, . . .N\ — 1}, 
sampled at (Sti)^ 1 = 4096Hz, consisting of Ni = 
T/8t\ = 245760 points. This is obtained by down- 
sampling the raw data stream by a factor of 4. 
The data are downsampled to 4096 Hz rather than 
2048 Hz to ensure that the rolloff of the associ- 
ated anti-aliasing filter is outside the frequency 
range being analyzed. The raw LI data are related 
to gravitational-wave strain by the calibration re- 
sponse function R\{f) constructed as described in 

sec.inng 

• For Al, a complex time series {s 2 [k]\k = 0, . . . N 2 — 
1}, sampled at (<5^2) 1 = 250 Hz, consisting of 
N 2 = T/St 2 = 15000 points. This is cal ibrated 
to represent strain data as described in Scc. lIIIB~2l 
but still heterodyned. 

To produce an approximation of the Fourier transform 
of the data from detector i, the data are multiplied by 
an appropriate windowing function, zero-padded to twice 
their original length, discrete-Fourier-transformed, and 
multiplied by SU. In addition, the LI data are multiplied 
by a calibration response function, while the Al data 
are interpreted as representing frequencies appropriate 
in light of the heterodyne. For LI, 



Sl(fl) ~ si 



JVi-l 



Mh) E «"ib"]9ib"]exp( 4^ ) 



3=0 



2N 1 

l=-N x ,...,Ni-\ (4.12) 



where fa = ^ is the frequency associated with the £th 
frequency bin. In the case of Al, the identification is 
offset by l\ = f h • (2T) = 107880: 

Mft)~«2[£\ 

rn rn ( -i2k {I - l h 2 )k\ _ + 
:= w 2 [k}q 2 [k}cxp ( ^ ±i- ) St 2 



k=0 



2N 2 



t%-N 2 ,...,t%+N 2 -l (4.13) 



As is shown in [361 ]. if we construct a cross-correlation 
statistic 



Y : = 



(4.14) 



the expected mean value in the presence of a stochastic 
background with spectrum S gw (f) is 



2. Discrete-Time Cross-Correlation 

Explicitly, the time-series inputs to the analysis 
pipeline, from each T = 60 sec of analyzed data, are: 



/i := (Y) « wm § E -jp [Q(m lift) Sgw(//) 

£=« mm 

(4-15) 

where W\W 2 is an average of the product of the two win- 
dows, calculated using the points which exist at both 
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sampling rates; specifically, if r\ and r-i are the smallest 
integers such that St 1 /St 2 = r 1 /r 2 = 125/2048, 



N/(nr 2 )-l 

nr 2 \ - 



w\w 2 = —rr- / Wi[nr2]w2[nri] 



N 



=o 



(4.16) 



Note that while the average value given by ()4.15j) is mani- 
festly real, any particular measurement of Y will be com- 
plex, because of the band-pass associated with the het- 
erodyning of the Al data. As shown in [36| . the real 
and imaginary parts of the cross-correlation statistic each 
have expected variance 



a 2 := \{Y*Y) 



T ■ 



l2 



qui) PidftDmm 



(4.17) 



where once again wfw 2 is an average over the time sam- 
ples the two windows have "in common" : 



2„,,2 



\ w 2 



nr 2 
N 



N/{ ri r 2 ) 



J2 (w l [nr2\f{w 2 [nr 1 ]) 2 (4.18) 



3. Construction of the Optimal Filter 



To perform the cross-correlation in (|4.14| . we need to 
construct an optimal filter by a discrete approximation 
to (|4.9p . We approximate the power spectra P\{f) using 
Welch's method [37} ; as a consequence of the averaging of 
pcriodograms constructed from shorter stretches of data, 
the power spectra are estimated with a frequency reso- 
lution Sf which is coarser than the 1/2T associated with 
the construction in (|4.14[) . As detailed in [22j, we han- 
dle this by first multiplying together [si(fe)]* and S2(fe) 
at the finer frequency resolution, then summing together 
sets of 2T Sf bins and multiplying them by the coarser- 
grained optimal filter. For our search, Sf = 0.25 Hz and 
T = 60 sec, so 2TSf = 30. 



4- Power Spectrum Estimation 

Because the noise power spectrum of the LLO can vary 
with time, we continuously update the optimal filter used 
in the cross-correlation measurement. However, using an 
optimal filter constructed from power spectra calculated 
from the same data to be analyzed leads to a bias in the 
cross-correlation statistic Y, as detailed in To avoid 
this, we analyze each T = 60 sec segment of data using an 
optimal filter constructed from the average of the power 
spectra from the segments before and after the segment 
to be analyzed. This method is known as "sliding power 
spectrum estimation" because, as we analyze successive 



segments of data, the segments used to calculated the 
PSDs for the optimal filter slide through the data to re- 
main adjacent. The Sf = 0.25 Hz resolution is obtained 
by calculating the power spectra using Welch's method 
with 29 overlapped 4-second sub-segments in each 60- 
second segment of data, for a total of 58 sub-segments. 



V. POST-PROCESSING TECHNIQUES 

A. Stationarity Cut 

The sliding power spectrum estimation method de- 
scribed in Sec. HVB~4"l can lead to inaccurate results if the 
noise level of one or both instruments varies widely over 
successive intervals. Most problematically, if the data 
are noisy only within a single analysis segment, consid- 
eration of the power spectrum constructed from the seg- 
ments before and after, which are not noisy, will cause 
the segment to be over-weighted when combining cross- 
correlation data from different segments. To avoid this, 
we calculate for each segment both the usual estimated 
standard deviation 07 using the "sliding" PSD estimator 
and the "naive" estimated standard deviation a'j using 
the data from the segment itself. If the ratio of these two 
is too far from unity, the segment is omitted from the 
cross-correlation analysis. The threshold used for this 
analysis was 20%. The amount of data excluded based 
on this cut was between 1% and 2% in each of the three 
orientations, and subsequent investigations show the final 
results would not change significantly for any reasonable 
choice of threshold. 



B. Bias-Correction of Estimated Errorbars 

Although use of the sliding power spectrum estima- 
tor removes any bias from the optimally-filtered cross- 
correlation measurement, the methods of [34[ show that 
there is still a slight underestimation of the estimated 
standard deviation associated with the finite number of 
periodograms averaged together in calculating the power 
spectrum. To correct for this, we have to scale up the er- 
rorbars by a factor of 1 + l/(A^ avg x 9/11), where N avg is 
the number sub-segments whose periodograms are aver- 
aged together in estimating the power spectrum for the 
optimal filter. For the data analyzed with the sliding 
power spectrum estimator, 29 overlapping four-second 
sub-segments are averaged from each of two 60-second 
segments, for a total A^ avg = 58. This gives a correction 
factor of 1+1/(58x9/11) = 1.021. The "nai've" estimated 
sigmas, derived from power spectra calculated using 29 
averages in a single 60-sccond segment, are scaled up by 
a factor of 1 + 1/(29 x 9/11) = 1.042. 
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C. Combination of Analysis Segment Results 

As shown in (l9j , the optimal way to combine a series 
of independent cross-correlation measurements {Yj} with 
associated one-sigma errorbars is 



v o P t _ Ej°7 y i 

~ _-2 



(Tyopt 



-1/2 



(5.1a) 
(5.1b) 



To minimize spectral leakage, we use Hann windows in 
our analysis segments, which would reduce the effective 
observing time by approximately one-half, so we over- 
lap the segments by 50% to make full use of as much 
data as possible. This introduces correlations between 
overlapping data segments, which modifies the optimal 
combination slightly, as detailed in [35| . 



D. Statistical Interpretation 

The end result of the analysis and post-processing of 
a set of data is a an optimally combined complex cross- 
correlation statistic Y° pt with a theoretical mean of R T 
and an associated standard deviation of ay°p* on both the 
real and imaginary parts. We can construct from this our 
overall point estimate on the unknown actual value of Cl R 
and the corresponding one-sigma crrorbar: 



Qr = Y opt /T 
CQ = cry opt /T 



(5.2a) 
(5.2b) 



For a given value of Ojj, and assuming <tq to be given, 
the likelihood function for the complex point estimate to 
have a particular value VL R = x + iy is 



P(x,y\Q R ,<rn) 



d 2 P 1 
dx dy 27T(7jj 

1 



\x + iy - Or| 
24 



(5.3) 



Given a prior probability density function on fl R , Bayes's 
theorem allows us to construct a posterior 



P(fl R \x,y,(j n ) 



= P{x,y\n R ,a n )P(n R ) 

P(x,y\<rn) (5.4) 
K e -(«-n«) a /aoSp(n H ) 



where x = KeQ R . In this work we choose a uniform 
prior over the interval [0,f2 max ], where O max is chosen 
to be 116 (the previous best upper limit at frequencies 
around 900 Hz [20j). except in the case of the hardware 
injections in Sec. IVII C[ where the a priori upper limit is 
taken to be well above the level of the injection. 



Given a posterior probability density function (PDF), 
the 90% confidence level Bayesian upper limit fiuL is 
defined by 



dn R P(Sl R \x,y,a n ) =0.90 



(5.5) 



Alternatively, any range containing 90% of the area under 
the posterior PDF can be thought of as a Bayesian 90% 
confidence level range. To allow consistent handling of 
results with and without simulated signals, we choose 
the narrowest range which represents 90% of the area 
under the posterior PDF. (This is the range whose PDF 
values are larger than all those outside the range.) For a 
low enough signal-to-noisc ratio HcQ r /<7q, this range is 
from to r^uL- 



E. Treatment of Calibration Uncertainty 

In reality, the conversion of raw data from the LIGO 
and ALLEGRO GW detectors into GW strain is not per- 
fect. We model this uncertainty in the calibration process 
as a time- and frequency-independent phase and magni- 
tude correction, so that fl R = x + iy and ctq are actually 
related to fl R e A+l ^, where A and <f> are unknown ampli- 
tude and phase corrections; the likelihood function then 
becomes 



P(x,y\tt R ,an,A, 



1 



2^ 



( \x + iy-n R e A+i 'f'\ 
6XP 



(5.6) 

Given a prior PDF P(A, <fi) for the calibration correc- 
tions, we can marginalize over these nuisance parameters 
and obtain a marginalized likelihood function 

P(x,y\n R ,an) 

dA J dct>P{x,yp R ,a u ,AA)P{K4>) (5.7) 

We take this prior PDF to be Gaussian in A and <f>, with a 
standard deviation added in quadrature from the quoted 
amplitude and phase uncertainty for the two instruments. 
For LI, this is 5% in amplitude and 2° in phas e |39l | and 
for Al it is 10% in amplitude and 3° in phase [211 ] . 



VI. ANALYSIS OF COINCIDENT DATA 

A. Determination of Analysis Parameters 

To avoid biasing our results, we set aside approxi- 
mately 9% of the data, spread throughout the run, as 
a playground on which to tune our analysis parameters. 
Based on playground investigations, we settled on the 
following parameters for our analysis: 

• Overlapping 60-second Hann-windowed analysis 
segments 
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• Frequency range 850 Hz < / < 950 Hz, 0.25 Hz res- 
olution for optimal filter 

• LI data downsamplcd to 4096 Hz before analysis 

• Removal of the following frequencies from the op- 
timal filter: 900 Hz (2.25 Hz wide), 904 Hz (0.25 Hz 
wide). 

The frequencies to remove were chosen on the basis of 
studies of the coherence of stretches of LI and Al data 
(see Fig. [2J; 900 Hz is the only harmonic of the 60 Hz 
power line within our analysis band, and detectable co- 
herence is seen for frequencies within 1 Hz of the line. 
904 Hz is notched out because, as the heterodyne fre- 
quency, it corresponds to DC in the heterodyned Al data. 
After completing the cross-correlation analysis, we com- 
puted the coherence from the full run of data, as shown 
in Fig. [SI The results are similar to those from the play- 
ground, except for a lower background level, and a feature 
at the heterodyning frequency. 

The relevant range of frequencies can be determined by 
looking at the support of the integrand of (|4.1ip . known 
as the sensitivity integrand. The overall sensitivity in- 
tegrand, constructed as a weighted average over all the 
non-playground data used in the analysis, is plotted in 
Fig. [6l The area under this curve for a range of frequen- 
cies is proportional to that frequency range's contribution 
to cr~ 2 . We see that the integrand does indeed become 
negligible by a frequency of 850 Hz on the lower end and 
950 Hz on the upper end. We further see that most of 
the sensitivity comes from a 20-Hz wide band centered 
around 915 Hz. 



B. Cross-Correlation and Upper Limit Results 

After data quality cuts, exclusion of the "playground" , 
and application of the stationarity cut described in 
Sec. IV Al 44806 one-minute segments of data were ana- 
lyzed, for an effective observing time of 384.1 hours (con- 
sidering the effects of Hann windowing), of which 181.2 
hours was in the XARM orientation, 114.7 in the YARM 
orientation, and 88.2 in the NULL orientation. The re- 
sults are shown in Table [Til No statistically significant 
correlation is seen in any orientation, and optimal combi- 
nation of all data leads to a point estimate of 0.31 + 0.30i 
for with a one-sigma errorbar of 0.48 each on the 
real and imaginary parts. 

The results in Table [UJ include the ORF describ- 
ing the geometry in the optimal filter. This means 
an orientation-independent non-GW cross-correlation 
present in the data would change sign between XARM 
and YARM, and would look much larger in the NULL 
result. One way to remove the effects of the observing 
geometry and compare orientation-independent cross- 
correlations is to remove the j(f) from (|4.9p . Since the 
ORF for each orientation is nearly constant across the 
observing band, and notably across the region of peak 
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FIG. 4: LLO-ALLEGRO (Ll-Al) coherence, calculated from 
48.66 hours of playground data spanning nearly 30 days. The 
only significant feature is the power line harmonic at 900 Hz. 
The closeup view in the second plot shows that the coherence 
is insignificant beyond 1 Hz away from the line. Based on this, 
we mask out the nine 0.25-Hz frequency bins around 900 Hz 
from our analysis. 



sensitivity, it is sufficient to multiply the overall results 
in each case by 7(915 Hz). This is shown in Table [TTTl 
where we again see no significant cross-correlation, and 
sensitivities whose relative sizes are well explained by the 
differing observing times. 

We use the methods of Sec. IVDI and IV El to construct 
a posterior PDF from the overall cross-correlation mea- 
surement of 0.31 + 0.30i and estimated errorbar of 0.48, 
taking into consideration the nominal calibration uncer- 
tainty of 11% in magnitude and 3.6° in phase to obtain 
the posterior PDF shown in Fig. \7\ The narrowest likely 
90% confidence interval on H,r is [0, 1.02]. We thus set an 
upper limit of 1.02 on flji = il gw (fji), which translates to 
an upper limit on ^(915 Hz) of (1.5 x 10" 23 Hz~ 1/2 ) 2 . 
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A1 L1 S4 full-run Coherence — combined orientations 
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FIG. 5: LLO- ALLEGRO (Ll-Al) coherence, calculated from 
all S4 data without regard to playground status. Again, the 
900 Hz line is seen to be comfined to a 2 Hz wide range. Ad- 
ditionally, a feature at the heterodyning frequency of 904 Hz 
(which was masked a priori in our main cross-correlation anal- 
ysis) becomes visible. 



VII. VALIDATION VIA SIGNAL INJECTION 

To check the effectiveness of our algorithm at detecting- 
stochastic GW signals, wc performed our search on data 
with simulated waveforms "injected" into them. This was 
done both by introducing the simulated signals into the 
analysis pipeline (software injections), and by physically 
driving both instruments in coincidence (hardware injec- 
tions). Hardware injections provide an end-to-end test of 
our detection pipeline and also a test on the calibration 
accuracies of our instruments, but are necessarily short in 
duration because they corrupt the GW data taken during 
the injection. Software injections can be carried out for 
longer times and therefore at lower signal strengths, and 
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FIG. 6: The sensitivity integrand for the data used in the 
cross-correlation analysis, normalized so its integral equals 
unity. The area under this curve, between two frequencies, is 
the fractional contribution to a~ 2 from that range of frequen- 
cies. Notice that the nine frequency bins masked out around 
900 Hz and the one at 904 Hz give no contribution, and that 
the sensitivity integrand is also suppressed at other frequen- 
cies corresponding to lines in the Al noise power spectrum. 



Type 


Teff 

(hrs) 


Point Estimate 


Error Bar 


XARM 


181.2 


0.61 + 0.25i 


0.56 


YARM 


114.7 


-0.47 + 0.47i 


0.90 


non-NULL 


295.8 


0.31 + 0.3H 


0.48 


NULL 


88.2 


10.96 - 43.89i 


28.62 


all 


384.1 


0.31 + 0.30i 


0.48 



TABLE II: Results of optimally-filtered cross-correlation of 
non-playground data. Results are shown for data in each of 
three orientations (XARM, YARM, and NULL). Addition- 
ally, the XARM and YARM results are combined with the 
optimal weighting (proportional to one over the square of the 
errorbar) to give a "non-NULL" result, and results from all 
three orientations are optimally combined to give an overall 
result. In each case, T e fj is the effective observing time includ- 
ing the effects of overlapping Hann windows. Note that since 
the non-NULL data are much more sensitive than the NULL 
data, they dominate the final result. Note also that because 
the optimal filter includes the the ORF, the relative orien- 
tation of LLO and ALLEGRO is already included in these 
results. This is reflected, for example, in the large errorbars 
on the NULL result. 



can be repeated to perform statistical studies. Software 
injections, however, cannot check for calibration errors. 
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Teg 






Type 


(hrs) 


Point Estimate 


Error Bar 


XARM 


181.2 


0.58 + 0.24i 


0.53 


YARM 


114.7 


0.42 - 0.42i 


0.80 


NULL 


88.2 


0.35 - 1.40« 


0.91 



TABLE III: The cross-correlation results of Table ILT1 scaled 
by 7(915 Hz) from Table [H the value of the ORF at the fre- 
quency of peak sensitivity. This gives a sense of the "raw" 
cross-correlation, independent of the orientation-dependent 
geometrical factor. The different observing times explain the 
remaining variation in the one-sigma errorbar for the measure- 
ment, which should be inversely proportional to the square 
root of the observing time. 



Posterior PDF & 90% conf band from all non-PG data 




FIG. 7: Posterior probability density function associated with 
the overall combined point estimate of 0.31 + 0.30? and esti- 
mated errorbar of 0.48, considering the uncertainty in the 
phase and amplitude of the calibration. The shaded region 
represents 90% of the area under the curve, leading to an up- 
per limit on f2_R of 1.02, which corresponds to a gravitational 



wave strain of 1.5 
915 Hz. 



x 10 iA Hz 1/2 at the peak frequency of 



A. Signal Simulation Algorithm 

To simulate a correlated SGWB signal in an interfer- 
ometer and a bar, the formulas in, e.g., [l9| need to be 
generalized slightly. This is because the ORF of a detec- 
tor with itself is in general [27j 



d ab d n 



1 



(7.1) 



which is unity for an IFO with perpendicular arms 
but 4/3 for a bar (|2.3p . Writing this quantity for detector 
1 or 2 as 7n or 722, respectively, makes the required 



cross-correlations in a simulated SGWB signal 

(K(f)hi(f)) = \Kf-f)S ew (f) in (7.2a) 
(K(f)Mf')) = \s(f-f)S gvJ (f) 112 (f) (7.2b) 
<W)M/')> = \s(f - /') Sgw(/) 722 . (7.2c) 

The above expressions do not determine a unique algo- 
rithm for converting a set of random data streams into 
individual detector strains. One possible prescription is 



hi(f) = 
Mf) = hi(f) 



s g w(/) V7irM/) + «M/)) 

712 (/) 



(7.3a) 



7n 



+ 2 \/ ^sw If) ( 722 



MR 

711 



(x 2 (f)+iy2(f)) , 
(7.3b) 



where x\{f), Di{f), £2(/), and yi (/) are statistically in- 
dependent real Gaussian random variables, each of zero 
mean and unit variance. In the above pair, 712(7) is 
used only in the construction of h<2(f) and not of hi(f). 
A different pair, where 712 is explicitly included in the 
calculation of both strains more symmetrically, can be 
defined as follows: Let z k (f) := (x k (f) + iy k (f))/V2 be 
a pair (fc = 1,2) of complex random functions and let 

s := \J\ —712/(711722)- Then, the second pair can be 
expressed as: 



hxU) = 



Sgw(/) 



/7TT(a(/)z 1 (/) + 6(/)z 2 (/)) (7.4a) 



(7.4b) 



where a = + s)/2 and b = 712/1/2(1 + s)7n 722 are 
determined completely by the three ORFs. Either pair 
of simulated strains obeys (|7.2p . The signals for software 
injections were generated using (|7.4[) ; those for hardware 
injections were generated by an older code which used 
(17. 3D. Further details of simulated signal generation arc 
in~|4|. 



B. Results of Software Simulation 

We performed software injections into the full S4 co- 
incident playground, 4316 overlapping one-minute anal- 
ysis segments with an effective observing time of 37.0 
hours considering the effects of Hann windows (16.7 hours 
of this is in the XARM orientation, 11.2 hours in the 
YARM orientation, and 9.0 hours in the NULL orienta- 
tion). We injected constant-f2 gw (/) spectra of strengths 
corresponding to fin —1.9, 3.9, 9.6, and 19, as well as a 



14 



Qr injected 


Point Estimate 


Error Bar 


90% conf int 





0.32 - l.OOi 


1.54 


[0.00,2.74] 


1.9 


2.22 - 0.86i 


1.55 


[0.09,4.35] 


3.9 


4.14-0.79; 


1.55 


[1.61,6.66] 


9.6 


9.89 - 0.65i 


1.56 


[7.32,12.45] 


19 


19.56 - 0.49i 


1.58 


[16.96,22.15] 



TABLE IV: Results of software injections. All figures are for 
constant- f2 gw (/) and listed by Hr level. The 90% confidence 
level ranges are calculated without marginalizing over any 
calibration uncertainty. 



test with the SGWB amplitude set to to reproduce 
the analysis of the playground itself. Note that even 
the strongest of these injections does not produce cor- 
relations detectable in an individual one-minute analysis 
segment. The results are summarized in table IIVI In 
each case, the actual injected value of Qr is consistent 
with the real part of the point estimate to within the 
one-sigma estimated error bar; the imaginary part of the 
point estimate remains zero to within the errorbar. The 
results for injections at Qr = 1.9 and stronger would cor- 
respond to statistical "detections" at the 90% or better 
confidence level. 



C. Results of Hardware Injection 

During S4, a set of simulated signals was injected in 
the hardware of ALLEGRO and LLO. These injections 
served to test the full detection pipeline as well as the 
calibrations of both instruments. As described in detail 
in [40j | , the preparation of simulated waveforms for hard- 
ware injections requires application of the transfer func- 
tion of the hardware component that is actuated, such as 
one of the two end test masses in LLO, to the theoreti- 
cal strain for that instrument. Subsequent refinements to 
instrumental calibration mean that the precise injected 
signal strength is determined after the fact. Six hardware 
injections performed during S4, each 1020 seconds long, 
had an effective constant f2 gw (/) of 8100. The series of 
injections we call "A" and "B" were performed during 
the XARM and NULL observation periods, respectively. 
Independent of the physical orientation of ALLEGRO, 
the injection and analysis was performed for three differ- 
ent assumed orientations, producing "plus" (aligned, as 
in the XARM orientation), "minus" (anti-aligned, as in 
YARM), and "null" (misaligned, as in NULL) injections 
in each series. The results (after correcting for known 
phase offsets in the injection systems) for each of the six 
injections are shown in Table IV1 

The results show some variation of magnitude and 
phase of the point estimates, especially for the "null" in- 
jections. However, all the injection results are consistent 
with the injected signal strength to within statistical and 
systematic uncertainties. This is illustrated informally 



in Fig. which shows the point estimates on the com- 
plex plane, each surrounded by an error circle of radius 
2.15 times the corresponding estimated errorbar. (This 
radius was chosen because 90% of the volume under a 
two-dimensional Gaussian falls within a circle of radius 
2.15(7.) Those circles all overlap with a region centered at 
the actual injection strength illustrating the magnitude 
and phase uncertainty in the calibration. The system- 
atic error can be more quantitatively evaluated using the 
method of Sec. IV El to produce a posterior PDF associ- 
ated with each injection measurement. This is illustrated 
for the optimally-combined point estimate of 7448 + 65?' 
and associated estimated errorbar 47 in Fig.[9j and ranges 
corresponding to the most likely 90% confidence range 
under the posterior PDFs (with and without marginal- 
ization) are included in Table [V] For each of the six 
injections, as well as for the combined result, the actual 
injected value of 8100 falls into the range after marginal- 
ization over the calibration uncertainty. 



VIII. COMPARISON TO OTHER 
EXPERIMENTS 

The previous most sensitive direct upper limit at the 
frequencies probed by this experiment was set by cross- 
correlating the outputs of the EXPLORER and NAU- 
TILUS resonant bar detectors (2(j. They found an up- 
per limit on /if 00 f^ gw (907.20 Hz) of 60. Using the value of 
72km/s/Mpc for the Hubble constant, that translates to 
a limit of 116 on fl gw (907.20 Hz), upon which our limit 
of 1.02 is a hundredfold improvement. 

Data from LLO, taken during S4, were also corre- 
lated with data from the LIGO Hanford Observatory 
(LHO) to set an upper limit on f2 gw (/) at frequencies be- 
tween 50 Hz and 150 Hz (24[. Correlations between LLO 
and LHO are not suited to measurements at high fre- 
quencies because of the effects of the ORF, illustrated 
in Fig. [2] For comparison, rough measurements using 
S4 LLO-LHO data and a band from 850 Hz to 950 Hz 
yield upper limits of around 20, while those confined to 
905 Hz < / < 925 Hz (the band contributing most of the 
Ll-Al sensitivity give upper limits of around 80. 

Correlations between the 4km and 2km IFOs at LHO, 
known as HI and H2, respectively, are not suppressed 
by the ORF, which is identically unity for colocated, 
coaligned IFOs. Since HI and LI have comparable sensi- 
tivities, the most significant factor in comparing H1-H2 
to Ll-Al sensitivity is the relative sensitivities of Al and 
H2. Since H2 was about a factor of 50 (in power) more 
sensitive than Al during S4, averaged across the band 
from 905 Hz to 925 Hz, we would expect an H1-H2 cor- 
relation measurement during S4 to be a factor of 7 or 
better more sensitive than Ll-Al as a measure of f2 gw (/) 
at these frequencies. However, the fact that HI and H2 
share the same physical environment at LHO necessitates 
a careful consideration of correlated noise which is ongo- 
ing @. 
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Injection 


7^fl 
Error Bar 


7^fl 
Point Estimate 


n R 

Error Bar 


£Ir Point 
Value 


Estima 
Mag 


te 

Phase 


unmj 
min 


irg. range 
max 


marg 
min 


. range 
max 


A- minus 


83 


-6623 - 126i 


93 


7403 + 141i 


7404 


1.1° 


7250 


7555 


6212 


8820 


A-null 


99 


205 + 19i 


3106 


6429 + 607i 


6457 


5.4° 


1565 


11294 


1435 


11562 


A-plus 


83 


6983 + 64i 


87 


7325 + 67i 


7325 


0.5° 


7182 


7468 


6146 


8726 


B-minus 


95 


-6845 + 49i 


106 


7650 — 55i 


7651 


—0.4° 


7478 


7825 


6417 


9115 


B-null 


111 


366 - 50i 


3486 


11492 - 1576i 


11600 


-7.8° 


5950 


17035 


5857 


17365 


B-plus 


92 


7128 + 77i 


96 


7477 + 80i 


7477 


0.6° 


7317 


7634 


6272 


8907 


all 


N/A 


47 


7448 + 65i 


7448 


0.5° 


7371 


7526 


6256 


8867 



TABLE V: Results of hardware injections. Simulated waveforms with an effective signal strength of fi gw (/) = 8100 were 
injected coincidentally in the ALLEGRO and the LIGO-Livingston (LLO) detectors. The "A" and "B" sets of injections took 
place during the XARM and NULL observing times, respectively. Independent of the actual orientation, the simulated signals 
were generated and analyzed assuming different orientations: YARM, NULL, and XARM orientations were assumed for the 
injections labelled "minus" , "null" , and "plus" , respectively. The first pair of columns shows the errorbars and point estimates 
scaled by 7(915 Hz) to give a "raw" cross-correlation as in Table|IIT] (Note that since the "null" alignment represents not-quite- 
perfect misalignment, as noted in Tab le|T] [7(915 Hz) = 0.03 rather than zero], the injection still leads to a statistically significant 
cross-correlation even in the "null" orientation.) Note that the errorbars, thus scaled, are comparable for all six injections, 
while the level of correlation or anti-correlation depends on the orientation associated with the injection being analyzed. The 
subsequent columns relate to the standard cross-correlation statistic, with the ORF included in the optimal filter, so the 
relative insensitivity in the "null" alignment is reflected in large errorbars, while the point estimates are all positive and in the 
vicinity of the injected value of 8100. All point estimates given include corrections for known phase offsets associated with the 
injection system. The row labelled "all" gives the optimally-weighted combination of all six results. The point estimates and 
one-sigma estimated errorbars were combined to give 90% confidence ranges with and without marginalization over calibration 
uncertainty, mimicking the statistical analysis described in Sec. IV Dl and fVEl 



Work is also currently underway to search for a SGWB 
at frequencies around 900 Hz by correlating data from 
the Virgo IFO with the resonant bar detectors AURIGA, 
EXPLORER, and NAUTILUS Q. 

Finally, an indirect limit can be set on SGWB strength 
due to the energy density in the associated gravitational 
waves themselves, which is given by 



Pgw — Pcrit 



^gw(/) 
/ 



df 



The most stringent limit is on a cosmological SGWB, set 
by the success of big-bang nucleosynthesis, is /c gw /Pcrit < 
1.1 x 10 -5 0. In comparison, a background of the 
strength constrained by our measurement, — 1.02, 
would contribute about 2 x 10~ 2 to /? gw /A;rit, if it were 
confined to the most sensitive region between 905 Hz 
and 925 Hz. (Spread over the full range of integration 
850 Hz < / < 950 Hz, it would contribute 1 x lO" 1 .) 
Note, however, that this nucleosynthesis bound is not 
relevant for a SGWB of astrophysical origin. 



roughly fivefold reduction in strain noise power at 900 Hz 
between S4 and S5, and because of the larger volume of 
data (roughly 20 times as much). Those two improve- 
ments could combine to lead to an improvement of about 
an order of magnitude in sensitivity. However, no 
immediate plans exist to carry out an analysis with S5 
data, because this incremental quantitative improvement 
in sensitivity would still leave us far from the level needed 
to detect a cosmological background consistent with the 
nucleosynthesis bound, or an astrophysical background 
arising from a realistic model. 

Additionally, the greater sensitivity of the H1-H2 pair 
means that a background detectable with Ll-Al would 
first be seen in H1-H2. In the event that a "surprise" cor- 
relation is seen in H1-H2 which cannot be attributed to 
noise, correlation measurements such as LLO- ALLEGRO 
and Virgo- AURIGA could be useful for confirming or rul- 
ing out a gravitational origin. 



X. CONCLUSIONS 



IX. FUTURE PROSPECTS 

LIGO's S5 science run began in November 5 with the 
aim of collecting one year of coincident data at LIGO 
design sensitivity. ALLEGRO has also been in opera- 
tion over that time period, so the measurement docu- 
mented in this paper could be repeated with S5 data. 
Such a measurement would be more sensitive due to Li's 



We have reported the results of the first truly het- 
erogeneous cross-correlation measurement to search for 
a stochastic gravitational-wave background. While the 
upper limit of 1.5 x 10 _23 Hz -1 ^ 2 on the strain of the 
SGWB corresponds to a hundredfold improvement over 
the previous direct upper limit on f2 gw (/) in this fre- 
quency band (2fj| . the amplitude of conceivable spectral 
shapes is already constrained more strongly by results at 
other frequencies [23|, [24[ . The lasting legacy of this work 
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FIG. 8: Visualization of hardware injection results. Each of the six point estimates of Qr is plotted on the complex plane, 
with an associated error circle of 2.15 times the estimated one-sigma errorbar. (This contains 90% of the volume under the 
corresponding likelihood surface.) The five-pointed star indicates the actual injected level of Qr = 8100. The dashed teardrop- 
shaped region indicates the calibration uncertainty, corresponding to a 2.15-sigma ellipse in log-amplitude/phase space. On 
the left we see that the two "null" injections are consistent in amplitude and phase with the actual injection, considering the 
statistical uncertainty associated with the real and imaginary parts of their point estimates. The plot on the right (in which 
the edge of the dashed calibration uncertainty teardrop can just be seen) shows that the "plus" and "minus" injections are 
all statistically consistent with each other, and consistent with the injection when systematic uncertainties associated with 
calibration are taken into account. 



is thus more likely the overcoming of technical challenges 
of cross-correlating data streams from instruments with 
significantly different characteristics. Most obviously, we 
performed a coherent analysis of data from resonant-mass 
and interferometric data, but additionally the data were 
sampled at different rates, ALLEGRO data were hetero- 
dyned and therefore complex in the time domain, and 
entirely different methods were used for the calibrations 
of both instruments. Lessons learned from this analysis 
will be valuable not only for possible collaborations be- 
tween future generations of detectors of different types, 
but also between interferometers operated by different 
collaborations. 
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